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Abstract 

A new expansion scheme to evaluate the eigenvalues of the generalized evolution operator 
(Frobenius-Perron operator) H q relevant to the fluctuation spectrum and poles of the 
order-g power spectrum is proposed. The "partition function" is computed in terms of 
unstable periodic orbits and then used in a finite pole approximation of the continued 
fraction expansion for the evolution operator. A solvable example is presented and the 
approximate and exact results are compared; good agreement is found. 
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1. Introduction 

The hallmark of chaos is the exponential growth of initial uncertainty, i.e., the sensi- 
tivity on initial condition. This stems from the fact that all periodic orbits are unstable, 
and the representative point wanders over the state space. As a result the dynamical be- 
havior of the system is quite irregular. This is the fundamental reason for the necessity of 
a statistical description of chaotic systems. 

It is interesting to know, however, that despite of their instability, many statistical 
quantities are well approximated by periodic orbits. The periodic-orbit approximation of 
the natural invariant measure has been introduced by Sinai, Bowen, Ruelle in the 1970's [1] 
in the context of hyperbolic systems. Cycles have then been used for the characterization 
of strange sets and local expansion rates [2] and for the calculation of damping rates in 
the conventional double time correlation function [3]. The main purpose of the present 
paper is to show how fluctuation and order-g power spectra, which have previously been 
introduced to describe the long-time statistics and the overall temporal correlations of 
chaotic fluctuations [4-7], can be calculated using periodic orbits and the continued-fraction 
expansion [6,8,9]. We will thus obtain a new expansion scheme for the eigenvalues of the 
generalized Frobenius- Perron operator H q , which plays an important role for the statistical 
description of chaotic fluctuations. 

The present paper is organized as follows. In section 2 we give a brief review of our 
approach to dynamical fluctuations. In section 3 we will derive fundamental equations 
which determine the eigenvalues in terms of unstable periodic orbits; they generalize the 
ones studied in Ref. 3. A new scheme for the estimation of the eigenvalues is proposed in 
section 4 by using the continued-fraction expansion. A solvable non-obvious example is 
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presented in section 5, and finite-pole approximations of the continued- fraction expansion 
are compared with the exact results. The paper is closed in section 6 with some final 
remarks. The relation between the generating function and the generalized evolution 
operator and the cycle expansion of the Fredholm determinant are discussed in appendices. 

2. Statistical characterization of dynamical fluctuations 

Consider a dissipative, chaotic <i-dimensional map 



Let Uj = u{x(j)} be an observed time series which is a unique, scalar function of x(j) 
generated by (2.1). The probability distribution p n (u) that the average 



x(j + l)=f(x(j)), 



X = (xi,X 2 , • ■ • ,X d ). 



(2.1) 




n-1 



(2.2) 



j=0 



takes a value u is asymptotically given by [4,6,10] 



p n (u) =< 5(u n - U) > 




(2.3) 



— s(u)n 



for large n. Here p*(x) is the natural invariant density, satisfying 




(2.4) 



with the Frobenius-Perron (FP) operator H. The generating function 



n-1 





asymptotically takes the form 

M q (n) ~ e^ )n , (2.6) 
for large n. The functions 4> and s are related to each other via [4,6,10] 

4>(q) = — mm[s(u) — qu]. (2.7) 

u 

The function <f>(q) plays a fundamental role for the overall characterization of temporal 
fluctuations. However it cannot describe the explicit temporal correlations in {u{x(j)}}. 
These can be singled out as follows. Let us define a weighted average 

u(q) = d(j){q)/dq 

(2.8) 

= lim («(x(0))e ? ™" (x(0)) )/M„(n) , 

n— »oo 

and the order-g power spectrum [7] 

I q (u) = lim (F n (u, x(0))d(u n (x(0)) - u{q))) / p n {u{q)) 

(2.9) 

= lim (F n (u;,a;(0))e 9n "™ (x W ) )/M g (n), 



where 



F n (uj,x) = - 
n 



n-l 

J2Hf J (x)}-u(q)]e-^ 

3=0 



(2.10) 



Then I q (oj) is identical to the power spectrum over time regions where the coarse-grained 
quantities u n take the same value u(q) and singles out the temporal correlation character- 
istic of the intensive parameter q. In addition, one can define the generalized double time 
correlation function C q (n) by an ensemble average as above. It turns out that I q (u) and 
C q (n) are related to each other via the Wiener-Khinchin theorem. 
If there exist only poles, one obtains [7] 

W = 'y><-> ^jf^ , (2.11) 



and the corresponding order-g double time correlation function takes the form 

C q (n) = £'tf<0 exp[-{ 7 « + iu q }n] , (2.12) 
I 

where {7^ } and {oj q } are respectively the sets of decay rates and characteristic frequen- 
cies. The ordinary power spectrum [11] is given by I q= o(ui) [7]. 

3. Determination of eigenvalues of H q by periodic cycles 

In appendix A (see also [5]), we show that eq. (2.5) can be rewritten as 

M q (n) = J H^p,(x)dx, (3.1) 

where H q is the generalized Frobenius-Perron (GFP) operator defined by [5,12] 

H q F(x) = j ' 8{x- f(y))e^F(y)dy ; (3.2) 

Hq is identical to the usual FP operator. Let v q and ipq be the eigenvalues and eigen- 
f unctions of H q , 

H q ^ q l \x) = ^\x). (3.3) 
Expanding p*(x) in terms of {ip q \x)}, and inserting it into (3.1) we obtain 

M q (n) = J24 l) H l) ] n ' ( 3 - 4 ) 
1 

with expansion coefficients {Jq }, (Yli^q = !)■ Here and throughout the paper we 
assume that H q has only discrete eigenvalues. Therefore, the characteristic function (f>(q) is 
determined by the largest eigenvalue u q ^ , which is positive because of the Perron theorem, 



viz. 



</>(?) = log *4 0) , (3.5) 



1 
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(tf> = 1). 

The order-g power spectrum and the double time correlation function can be expressed 
as in (2.11), the poles being determined by 

exp[ 7 «+z^)] = ^°)/^), (Z^O). (3.6) 

The coefficients in (2.11) are obtained from [7] 

K® = a q (0,l)a q (l,0) 

u{x}^(x) = Y,a q (ij')4 i 'n^, (3 ' 7) 

v 

assuming completeness of the eigenfunctions. The functions I q (uj) and C q (n) concern the 
explicit temporal correlation as singled out with the parameter value q from infinitely 
many aspects of temporal-correlation characteristics embedded in {u{x(j)}}. The ordi- 
nary double time correlation function C(n) = (5u{x (n)}5u{x (0)}) , (where 5u = u — (u)), 
simply describes only one form, corresponding to the limit q — > 0, (C(n) = Co(n)). 
Let us define the xy element of the operator H q by 

{ Hq )^ = 8{x-f{y))e^\ (3.8) 

and therefore 

n-1 

(^) xy = 6(x - r(y)) expfe u if J (y)}] ■ (3-9) 

3=0 

The partition function 



Z q (n) =TtH^= I (#£) xx dx , (3.10) 

is expanded as 



n-1 

'„" ' ' ]]] 7 x 

3=0 

exp[QE;=o«{/ J > (n) ' fc )}] 



/( i. _1_ 
6(x-f n (x))e X p[qY, U {f1(x)}]d 
a n 



(3.11) 



J 3 

k 



det[l - D n (x( n ^ k ) 
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Here the summation J2 k extends over all fixed points {x^ ,k }, where x^ ,k denotes the 
k-th fixed point of the n-th iterate f n {x), (k = 1, 2, • • • , n), 

x (n),k = fn( x (n),ky (3^) 

Furthermore (D n (x)) jk = d(f n (x))j/dx k , i.e., 

D n (x) = l[ + g(P(x)), (g( x )) kl = ^l, (3.13) 

Y[ + being the ordered product. 

The spectral decomposition of H q on the other hand yields 

i 

Together with (3.11) we thus obtain a set of fundamental relations between the eigenvalues 
{vq^} and the periodic cycles as contained in Z q (n), viz. 

5><'>] B = , (3.15) 
l 

(n = 1, 2, 3, • • •). The next problem is to solve for the {vq^} in terms of {Z q (n)} which are 
determined by periodic orbits as in (3.11). In Ref. [3] and appendix B, the cycle expansion 
method is used for this purpose. Our main aim here is to study the applicability of a 
continued fraction expansion, to which we turn now. 

4. Continued-fraction expansion for the eigenvalues 

Let A[/j] denote the Laplace transform of a function A(n), 

oo 

AMe^W 1 . (4.1) 



Applied to eq.(3.15) one finds 

^m = E(4t _/4) " 1 - (42) 

I U Q 

The eigenvalues {vq^} are thus determined by the poles of Z q \jj], which we estimate from 
a continued fraction expansion. 
For convenience, let us put 

(n)i = Z q (n), (4.3) 
(n = 1, 2, 3, • • •). Let (n)i obey the "equation of motion", 

n-l 

(n + l)i = 6i(n)i + J2( n ~ lhU)u (4-4) 

where b\ = (2)i/(l)i. In the memory term a new function (71)2 appears. By constructing 
an equation of motion for (71)2 similarly to (4.4), a new function (71)3 enters. Repeating 
this procedure successively, we have a set of functions {{n) k , k = 1, 2, • • •}, related to each 
other via 

n-l 

(n + l) fc = 6 fc (n) fe + J^(n - j)k+i(j)k, (4.5) 

where 6fc = (2) fc /(l) fc , (fc = 1, 2, • • •). Equation (4.5) is regarded as the equation of motion 
for (n)fc determined by the memory kernel (n)k+i- 
The Laplace transform of (n) k is given by 



n=l 

(1)* 



(4.6) 



1 - 6 fc /j - [/j] fc+ i/j 2 ' 
Thus we get the continued-fraction expansion [6] 

i-M 

1 - 62/x 2?T\ 

1 — 04fi — • • • 



Since {v q '} are identical to the reciprocal poles of Z q [fj], we find 



= V/4°, (4-8) 



where are poles of (4.7), i.e., the zeros of the denominator. Obviously, the continued- 

fraction representation (4.7) can be expanded as (4.2). 

Now let us turn to how to determine the coefficients {(l)fc, bk, k = 2, 3, 4, • • •} in (4.7). 
One should note that the quantity (n)k+i is determined from {(j)k,j = 1, 2, ■ • • , n + 2} by 



n-l 



(n) fc+ i = (TJ^U™ + 2 )fc - &fcfa + l)fc - ^(/)fc+i(n + 1 - /)*}, (4.9) 



with 

(l) fc+1 = J_{(3),-6 fc (2) fc }, (4.10) 

(fc = 1, 2, • • •). So the quantity (n)^ is given in terms of {(j)i, j = 1, 2, • • • , n + 2(k — 1)}. 
Especially, parameters (l)fc and (2)fc, (/c = 2, 3,4, •••), in (4.7) are given, respectively, by 
sets {Z q (j),j = 1,2, ■■ -,2k - 1} and = 1, 2, • • • , 2k}. How (l) fc and (2) fc are 

determined by {Z q (j)} is given by solving the recursion relation (4.9) with the initial 
condition (4.10). 

Since from a practical viewpoint the infinite continued-fraction expansion is intrac- 
table, we must develop a discard approximation. This can be carried out with the finite- 
pole approximation of the continued-fraction expansion as follows. Let us assume that the 
eigenvalues are ordered, is q °^ > l^q 1 ^ > l^ 2 '*! > " •> an d let us terminate (3.15) after m 
terms, i.e., 

m— 1 

£ H l) ] n = Z q (n). (4.11) 



This is reasonable, if the eigenvalues with / < m decay sufficiently rapidly. This is equiva- 
lent to the m-pole approximation (l) m+ i = 0,i.e., 

Mm = , (4-12) 

1 - b m [i 

and leads to 

M* = ^m+2-k m m > (fe=l,2,...,m). (4.13) 

Here coefficients are 

a fc (l) = l, (fe=l,2,...,m+l), (4.14) 

m 

a fc (2) = -^6„ (fe=l,2,...,m), (4.15) 

a m _i(3) = -6 m _ia m (2) - (l) m . (4.16) 

Furthermore 

a fc (Z) = a fc+ i(/) - b k a k+ i(l - 1) - (l)fc+iafc +2 (Z - 2), (4.17) 
for fe = 1, 2, • • • , m — 2 and 3 ^ m — k + 1, and 

ajfc(0 = -6 fc a fc+ i(/ - 1) - (l) fc+ ia fc+2 (/ - 2), (4.18) 

for I = m — k + 2 where fe = 1, 2, • • • , m — 2. 

The poles with the m-pole approximation are therefore the roots of the order-m 

algebraic equation 

Q m (fi)= J2 a i( l )v l ~ 1 = - ( 4 - 19 ) 
i=i 

By putting v = the equation Cl m (l/v) = is therefore identical to the characteristic 
equation of (3.3) with the m-eigenstate approximation. 

Let us comment on the validity of the finite-pole approximation of the continued- 
fraction expansion. It is obvious that the expansion (4.5) is meaningful only when {|6fc|} 



are not too large. Assume that for one /c-value, say k = K( <^m), {1)k vanishes (or, in 
a practical sense, is extremely small) for some g-value, q = Q. The quantity \bx\ then 
diverges at q = Q, if (2)k remains finite. Since a part of O m (//) contains terms h((x)bKfJ> 
at q = Q, h being a non-singular factor, the algebraic equation O m ( / u) = has a solution 
/U = at q = Q, i.e one reciprocal eigenvalue v diverges. As q is changed from below Q, 
one reciprocal eigenvalue gradually, say, increases, and diverges at q = Q, i.e. v = oo for 
q = Q — 0. For q = Q + 0, v = —oo. Near such points the continued-fraction expansion 
loses its validity. 

5. Example — Exact eigenvalues and finite-pole approximation — 

In this section we study a simple non-obvious solvable example and compare the finite- 
pole approximation of the continued-fraction expansion with exact results. The model we 
discuss is the Markov map 



where sq = a 1 and s± = a/(l — a). 

A. Partition function and eigenvalues 

First, we calculate the partition function Z q (n) from periodic orbits. Let Iq = [0, a] 
and I\ = (a, 1] be the dynamical partition of the full interval I = [0, 1], i.e., the subintervals 
on which / has constant slope. / maps Iq onto I and I\ onto Iq. Higher iterates of / 
give rise to finer subdivisions of I which can be labelled by strings of letters A and B, 
depending on whether the iterates come to lie in the subintervals Iq or ii respectively. 
Since I\ is mapped onto Iq, there can be no repetitions of the letter B, i.e., every B is 






(5.1) 



followed by an A. If the last letter of a string of length k is A, then the k-th iterate of the 
subinterval is mapped onto the full interval and thus has a fixed point (cycle of period k). 
If the letter is a 5, then there will be a fixed point only if x e Iq. 



of intermediate points in the intervals Iq and Ii, respectively. Since they depend only on 
the number and not on the ordering, periodic points with strings differing only in the 
ordering of symbols will give rise to the same derivatives and the same weights in the 
traces, eq. (3.11). 

We now turn to the computation of the number of points with the same derivatives. 
To this end, we consider the sum of slopes, Wx(k) : over all fixed points of f k (x). A fixed 
point of f k (x) on a subinterval is classified into two groups A and B according to the last 
letter of the string. Let W A {k), (a = or 1) denote the sum of slopes over the fixed points 
in the group A which are on I a . Similarly, W B (k) is defined for the fixed points in the 
group B. Then, we can write the recursion equation 



The derivatives of at its fixed points are Sq A s 



n B 
1 > 



where ua and ub are the number 



W A (k) = s [W A (k - 1) + W B (k - 1)] 1 



W B (k) = W A (k) = Sl W A (k-l) 



> 



(5.2a) 



W A (l) = s , W A (0) = 0. 



The sum Wx(k) can be written as 



W T (k) = W A (k) + W B (k) + W A (k) . 



(5.26) 



The recursion equation (5.2a) leads to 



W A {k) 



y+-y- 1 



y k +-y- , 



(5.3) 



where we put 



Inserting (5.3) and (5.4) into (5.26), we have 
W T (k)=y k + + y k _ 

i=0 v 7 j=0 w 

where [•] is Gauss's symbol. Because the slopes of f h (x) are written in the factorized forms 
sjsg - "' (j = 0, 1, 2, ... , [k/2]), the number of the fixed points with the same slope is given 
by the corresponding coefficient of the homogeneous equation (5.5). 

If we now specify the observable u{x} to be the characteristic function on the interval 
[a, 1], i.e., 

(0 ^ x ^ a) 

u{x} = 1 (5.6) 
' 1 s 1 (a<xZl), 



z » = Q " E ( 2 "J E (J) M'/W<ir J - 1)- (5.7) 



we can write the partition function (3.11) as 

, [»/2] 

" E 

i=0 v " 7 j'=0 

Now we evaluate eigenvalues of the operator if g with observable (5.6). As we are 
interested in an absolutely continuous measure, we consider the functional space T which 
consists of all linear combinations of polynomials on I Q and I\. Let x a , (a = or 1) denote 
the characteristic function of I a : Xcr(x) = 1 for x e I a , and for x I a . Then, a linear 
combination G GJF is defined by 

G(x) = Pj(x) X o(x) + P k (x)xi(x), (5.8) 

where Pj and Pk are polynomials on I. Because the map (5.1) is piecewise linear and the 
u{x} piecewise constant, the map H q :T— *T , 

(H q G)(x) = Is^Pjis^x) + s^Pk^x + a)e q ]xo(x) + Su 1 P J (so 1 x)xi(x) , (5.9) 



preserves the maximal degree of polynominals. For n = 0, 1, 2, . . ., the set of linear com- 
binations of polynomials of degree k (0 ^ k ^ n) is a subspace of J 7 , denoted by . 
The restriction of H q onto JF 1 ^ determines an operator -.T^ —tj 7 ^. As a linear 

combination G{x) in can be represented by a 2(n + l)-dimensional vector [11], and 
Hq n ^ can be written as a 2(n + 1) x 2(n + 1) matrix. Indeed, we have the upper-triangular 
matrix 



//4 0) 



q 



h 



(i) 



V 







(n) 

where each element of Hq is a 2 x 2 matrix and 



Si e 



\ 



(j = l,2,...,n+l). 



(5.10) 



(5.11) 



Eigenvalues of (1 ^ j ^ n + 1) are also eigenvalues of Hq 1 ' , and are given by 



(n) 



V ±V s o 2j +4e9(s si)-J 



(5.12) 



Since the full generalized Frobenius-Perron operator H q is obtained in the limit n — > oo, 
the eigenvalues of if g are 



,(0 



(j) 



+ 
(j) 



for / = 2j = 0, 2, 4, . . . 

for 1 = 2j + 1 = 1,3,5, .... 
Substituting (5.12) into (3.14) and summing up over / yields 



(5.13) 



Tri?? 

q \2 



E(;j(^)7W«s- J -i)- 



(5.14) 



i=0 ^ j=0 

We remark, that for this example the generating function (2.5) can be expressed in terms 
of eigenvalues of H^ only as 

M q (n) = {[1- ^} 2 H°r - [l - ^fMT}/^ " [2 " < ] ~ (5-15) 



B. Truncation of the continued- fraction expansion 

We will study the approximate methods to obtain eigenvalues of the GFP operator 
for the model (5.1) by truncating the continued-fraction expansion. Hereafter we consider 
three approximations with 4-, 6-, and 8-poles, for which we need traces up to order n = 8, 
10 and 12, respectively. The parameter is chosen as a = 0.75. 

Figure 1(a) shows the first twenty of the exact eigenvalues (5.13). Figures l(b-d) show 
the comparison between results of the finite-pole approximations and exact results. In each 
approximation there exist apparent deviations of approximate values from the exact results 
in certain regions of q. Such deviations are located near q- values where one or several 
coefficients (1)k vanish, thus leading to a divergence of \bk\ and a failure of the finite 
pole approximation of the continued fraction expansion. For the 4-pole approximation 
the region showing apparent deviation is rather wide. However, as the number of poles 
is increased, the number of such singular points increases, but the total area of effected 
(/-regions diminishes. In this way the approximate results approach the exact ones as the 
number of poles is increased. 

6. Concluding remarks 

In closing this paper let us add several remarks. The first point concerns the relation 
between the characteristic function <fi{q), especially for local expansion rates [13], and 
the decay rates {7^ } in I q {u>). Whereas 4>(q) describes the probability distribution of 
fluctuations, the decay rates {7^} characterize the explicit temporal correlations. (j)(q) 
is determined by the largest eigenvalue of H q and {7^ } by others [7]. In this sense, we 
conclude that decay rates generically have no relation with the characteristic function <p(q) , 



except in special cases. 

We proposed a continued-fraction expansion for the eigenvalues of the generalized 
Frobenius- Perron operator H q , whose expansion coefficients are determined in terms of 
periodic orbits of the dynamical equation. A solvable illustrative example was presented, 
and their eigenvalues were compared with those of truncated continued-fraction expansions. 
It should be noted that whereas the largest eigenvalue is insensitive to the approximation 
(except in several narrow regions of q, related to the divergence of the coefficients 
other eigenvalues which are relevant to the temporal correlations depent rather sensitively 
on the order of truncation in the continued fraction expansion. Except for this, the fi- 
nite pole approximation of the continued-fraction expansion turned out to be a powerful 
approximation. 

The ordinary double time correlation function is determined by the eigenvalues of the 
FP operator H [14]. The equations corresponding to (3.15) are 

£V? } ] n = Mn) ~ 1, (6.1) 
l 

where Yl'i stands for the summation except / = 0, because as a result of the existence of 
the invariant density the largest eigenvalue of {v^} is unity. As can be easily verified in 
our example, the eigenvalue vf 1 ' = t+ (eq. (5.13)) indeed passes through 1 for q — > 0. 
In Ref. 3, eq (6.1) was solved approximately using the cycle expansion of Ruelle zeta 
functions and Fredholm determinants dp- Because of possible poles in the Ruelle zeta 
functions, one has to use the full Fredholm determinant to determine the eigenvalues. As 
the calculations in appendix B show, in our solvable example, the Ruelle zeta functions 
are entire and the Fredholm determinant factorizes into a product of Ruelle zeta functions; 
thus there are no poles and the spectrum as computed from Ruelle zeta functions and the 
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Fredholm determinant coincides. This exact factorization of course is not noticable from a 
truncated cycle expansion of the Fredholm determinant as used in numerical calculations. 
The numerically determined eigenvalues behave very much as in the case of the continued 
fraction expansion, except that they show no poles. 

One might be lead to conclude from this example that the cycle expansion of appendix 
B seems to be better, but one can also find examples where the cycle expansion method 
has problems but where the continued fraction expansion works fine, e.g. in the case of the 
fully developped parabola. Thus the advantages and disadvantages of the two methods 
are more subtle, and certainly require further work. 

As reported in [7,8], {vq } can be derived as the poles of the Laplace transform of 
M q (n), where we need the explicit expression of moments M q (l), M q (2), M 9 (3), • • •, using 
the explicit information on the natural invariant density itself. The present approach 
however does not need the natural invariant density to obtain {vq }. 

Finally let us comment on how the order of truncation m should be chosen. Let r 
and 7 be the characteristic periodicity and the damping rate, respectively. Depending on 
the type of chaotic transitions and on the observables, r and 7 sometimes exhibit critical 
slowing-down. For example, for x^j), one component of the orbit x(t) itself, the damping 
rate 7 tends to zero with finite r as the control parameter approaches the band splitting 
point. And we find 7^0 and r — > 00 near the type I intermittency. Although there are 
two characteristic time-scales r and 7 _1 , it is expected that m should be chosen such that 
m ^> r, since r is relevant to periodicity. (See Ref.[15].) 
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Appendix A — Derivation of (3.1) — 

Equation (3.1) can be obtained by noting that 



p*(x) exp[<? ^2 u{f J (x)}]dx 
= / p,(x)e^Lexp[q^2u{f^x)}]dx, 



i=o 

where L is the time evolution operator, LG(x) = G(f(x)). Since f G\{x)LG2{x)dx = 
f[HG\(x)]G2(x)dx, eq. (A.l) can be written as 

/n-2 
[H q p*{x)\ exp[q ^ u{f{x)}]dx, (A.2) 

3=0 

where H q is the generalized Frobenius-Perron operator defined in (3.2). 
By repeating the above procedure successively, M q (n) is written as 

n— 3 



M q (n) = I [H 2 q p*(x)]exp[qJ2<f j (x)}}dx 

3=0 



(A.3) 



= J[H™ 1 p*(x)]exp[qu{x}]dx. 
By noting J G(x)dx = J HG{x)dx, eq.(A.3) leads to eq.(3.1). 

Appendix B — Cycle expansion and dynamical zeta function — 

The cycle expansion approach to solving the moment problem (3.15) consists again 
in forming the Laplace transform (4.1) so as to obtain the eigenvalues as inverses of the 
positions of poles. However, this Laplace transform is then expressed as a logarithmic 
derivative of a Zeta function, so that the eigenvalues are determined by zeros of this Zeta 
function. 



Consider again eq. (3.11), expressing the trace Z q (n) in terms of periodic orbits. Let 
p label all primitive cycles of period n p (n p is the smallest number of iterations after which 
the orbit closes). Note that both 

n-l 

u p = J2 u {f j ( x(n) ' k )} (B.l) 
j'=i 

and 

A p = \D n (x^ k )\ (B.2) 

are independent of the point where they are computed. A cycle of period n p has n p different 
periodic points. In the sum (3.11), it contributes whenever k = n p or a multiple thereof. 
Thus when summing on all periods in the Laplace transform (4.11), viz. 

oo 

Z q [^] = J2z q (n)^ n - 1 (B.3) 

71=1 

one can also express this as a sum on all primitive periodic orbits and their multiple 
traversals, 

zM^zi^-w^t^'- 1 - (BA) 

p r=l ' P' 

The usual manipulations (expansion of the denominator in a geometric series and summa- 
tion on r) [3] result in 

where the Zeta function 

oo oo 

= n n i 1 - ^ v exp^iApi- 1 v) = n c- 1 , (b&) 

V 3=0 j=0 

factorizes into dynamical zeta functions C -1 . As we will see, for the particular map stud- 
ied in section 5, each dynamical zeta function will be a quadratic polynominal and thus 
contribute two eigenvalues for each j. 
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The key to the evaluation of the products on periodic orbits is the topological organi- 
zation of the periodic orbits. In the example of section 5, there can be no two consequtive 
B's, since every point in the interval [a, 1] is mapped back into the interval [0, a]. This 
is the only exclusion rule. Defminig a new alphabet = A and 1 = BA, one finds that 
all strings of and 1 can be realized. Since moreover the maps is linear, all weights A p 
factorize according to 

K = A^A^b (B.7) 

where n and ri\ are the numbers of 0's and l's in the symbol string of p, respectively. The 
instabilities A p are 

= s , A AB = s si , (5.8a) 

and the observable (5.6) becomes 

u p = ri\ . (B.8b) 
One can then show [3,16] that the dynamical zeta functions become 

C" 1 = 1 -/"M" 1 ^ -^exp^lsosir^sosi)^'. (B.9) 

which is precisely the characteristic polynominal for the submatrices (5.11), C" 1 = det(l — 
^hq) for j = 0, 1, • • •. Thus in this case the cycle expansion yields the exact eigenvalues. 



- 20 - 



References 

[1] See, e.g., 

Sinai, Ya. G.: Uspekhi Math. Nauk 27, 27 (1972) 

Bowen, R.: Equilibrium States and the Ergodic Theory of Anosov Diffeomorphisms, 
Lecture Notes in Mathematics 470 (New York: Springer- Verlag, 1975) 
Ruelle, D.: Thermodynamic Fromalism in Encyclopedia of Mathematics and Its Ap- 
plications, Vol.5, (Addison-Wesley, Reading, 1978) 

[2] See, e.g., the following articles: 

Kai, T., Tomita, K.: Prog. Theor. Phys. 64, 1532 (1980) 
Takahashi, Y., Oono, Y.: Prog. Theor. Phys. 71, 851 (1984) 

Auerbach, D., Cvitanovic, P., Eckmann, J. -P., Gunaratne, G., Procaccia, I.: Phys. 
Rev. Lett. 58, 2387 (1987) 

Morita, T., Hata, H., Mori, H., Horita, T., Tomita, K.: Prog. Theor. Phys. 296 
(1988) 

Cvitanovic, P.: Phys. Rev. Lett. 61, 2729 (1988) 

Cvitanovic, P., Gunaratne, G. H., Procaccia, I.: Phys. Rev. A 38, 1503 
(1988) 

[3] Cvitanovic, P., Eckhardt, B.: J. Phys. A: Math. Gen. 24, L237 (1991) 
Christiansen, F., Paladin, G., Rugh, H. H.: Phys. Rev. Lett. 65, 2087 
(1990) 

Rugh, H. H.: Nonlinearity 5, 1237 (1992) 
Eckhardt, B.: Acta Phys. Pol. B 24, (1993) at press 

[4] Fujisaka, H., Inoue, M.: Prog. Theor. Phys. 77, 1334 (1987) 

- 21 - 



[5] Fujisaka, H., Inoue, M.: Prog. Theor. Phys. 78, 268 (1987) 

[6] Fujisaka, H.: in From Phase Transitions to Chaos (Topics in Modern Statistical 
Physics), eds. Gyorgyi, G., Kondor, I., Sasvari, L., and Tel, T., (World 
Scientific, 1992) 

[7] Fujisaka, H., Shibata, H.: Prog. Theor. Phys. 85, 187 (1991) 

Just, W., Fujisaka, H.: Physica D (1993), at press 
[8] Fujisaka, H., Inoue, M.: Prog. Theor.Phys. 78, 1203 (1987) 
[9] Just, W.: J. Stat. Phys. 67, 271 (1992) 

[10] Ellis, R. S.: Entropy, Large Deviations, and Statistical Mechanics, 
(New York: Springer- Verlag, 1985) 

[11] See, e.g., the following articles: 

Grossmann, S., Thomae, S.: Z. Naturforsch. 32a, 1353 (1977) 
Mori, H., So, B. C., Ose, T.: Prog. Theor. Phys. 66, 1266 (1981) 
Shigematsu, H., Mori, H., Yoshida, T., Okamoto, H.: J. Stat. Phys. 30, 
649 (1983) 

[12] Szepfalusy, P., Tel, T.: Phys. Rev. A 34, 387 (1986) 

[13] Badii, R., Heinzelmann, K., Meier, P. F., Politi, A.: Phys. Rev. A 37, 
1323 (1988) 

[14] Mori, H., et al. and Shigematsu, et al. in Ref. [11] 

[15] Just, W., et al. in Ref. [7] 

Kobayashi, T., Fujisaka, H., Just, W.: Phys. Rev. E, (1993) at press 

Just, W., et al, to be submitted 
[16] Artuso, R., Aurell, E., Cvitanovic, P.: Nonlinearity 3, 325 and 361 (1990) 

- 22 - 



Cvitanovic, P., Eckhardt, B.: Nonlinearity, (1993), at press 



- 23 - 



Figure Caption 

Fig.l (a) Exact eigenvalues of H q for the model (5.1) with a = 0.75. Twenty 
eigenvalues vs q are shown in the order of increasing values. The results from 
finite-pole approximations of the continued-fraction expansion are shown in 
figures (b) using 4 poles, (c) using 6 poles and (d) using 8 poles. Solid and 
dotted lines respectively indicate results of finite-pole approximation and exact 
results. As the number of poles is increased, the results tend to agree with the 
exact results. 
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